Phase diagram of a random-anisotropy mixed-spin Ising model 
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We investigate the phase diagram of a mixed spin- i -spin- 1 Ising system in the presence of 
quenched disordered anisotropy. We carry out a mean-field and a standard self-consistent Bethe- 
Peierls calculation. Depending on the amount of disorder, there appear novel transition lines and 
multicritical points. Also, we report some connections with a percolation problem and an exact 
result in one dimension. 
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I. INTRODUCTION 

Apart from their relevance to the description of fer- 
rimagnetic materials, mixed-spin models are interesting 
from a purely theoretical point of view, being among the 
simplest models to exhibit tricritical behavior. So they 
are especially convenient for studying the effects of in- 
homogeneities on the phase diagrams and the multidjit-, 
ical behavior of magnetic-systems. From a few exaclw'cl 
and several approximatemj calculations, we now have 
a good picture of the phase diagrams of mixed spin-i- 
spin-1 Ising models in the presence of a crystal field. Our 
aim in this work is to use this model Hamiltonian to in- 
vestigate the effects of disorder on the location of the 
transition lines and the tricritical point. 

The mixed-spin Ising model is described as a two- 
sublattice system, with spin variables cr = ±1 and 
S — 0, ±1, on the sites of sublattices A and B, respec- 
tively. Restricting the interactions to nearest-neighbors 
(belonging to different sublattices) and single-ion terms, 
the most general spin Hamiltonian in even spin space is 
written as 



H 



-J ^ (TiSj 



{ieA,jeB) 



(1) 



where the first sum is over nearest-neighbor pairs, the 
second sum is over the sites of sublattice B, and the pa- 
rameter J is assumed to be positive (ferromagnetic ex- 
change). For D > 0, the crystal field favors the Sj — 
states; the competition between anisotropy and ferro- 
magnetic exchange terms leads to the appearance of a 
tricritical point. It should be pointed out that one needs 
three parameters to describe the even space of the bet- 
ter known spin-1 Blume-Emery-Griffiths (BEG) model. 
In the present case of a mixed-spin Ising system, from 
the point of view of the calculations, the reduction of 
the analysis to a two-parameter space is a particularly 
attractive feature. 

There are exact calculations for the thermodynamic 
functions associated with the model Hamiltonian given 
by Eq. (pi) on a simple chain and on some three- 



coordinated two-dimensional structures. On a honey- 
comb lattice, the problem can be mapped onto a spin- 
i Ising model on a triangular lattice, which does not 
display a tricritical point .u This mixed-scin model can 
also be exactly solved on a Bethe latticeB (the deep in- 
terior of Cayley tree), leading to the^same results of a 
recent cluster-variational calculation.Q The results on a 
Bethe lattice with coordination q indicate the absence of 
a tricritical point for q < 5, as confirmed bju Migdal- 
Kadanoff renormalization-group calculations.Q In the 
infinite-coordination limit of the Bethe lattice, one re- 
gains the well-known results for the tricritical point dis- 
played by the Curie-Weiss (mean-field) version of thg 
model. An earlier approximate effective-field calculationQ 
predicted a tricritical poiaLfor q> A, but this result has 
been recently challengedJatl 

In order to analyze the effects of disorder, we consider 
the Hamiltonian 



H 



-J 



E 

{ieA,jeB) 



(2) 



where {Dj} is a set of independent, identically dis- 
tributed random variables associated with the binary 
probability distribution 



p{D,)=p5{D,) + {l-p)8{D,-D). 



(3) 



With this choice of disorder, and for D > qj, the ground 
state can be mapped onto a percolation problem in which 
the dilution affects the sites belonging to only one of the 
sublattices (corresponding to spin S* = 1). This associa- 
tion is easy to see if we note that a uniform crystal field 
D > qJ leads to Sj = for all j, which breaks the con- 
nectivity between the spin- 1/2 variables. The presence 
of a randomly located distribution of _D = crystal fields 
recovers that connectivity and, for sufficiently high val- 
ues of p, leads to the formation of a percolating cluster. 
In the rather artificial case of annealed disorder, cp,the 
honeycomb lattice, there is also an exact solutionE2l for 
the thermodynamic properties of the mixed-spin model 
described by Eqs. (H) and (^). It is interesting to remark 
that this solution in the annealed case reproduces the 
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critical concentration of the percolation problem associ- 
ated with the ground state of the model with quenched 
(frozen) disorder, which is equivalent to the usual perco- 
lation problem on the triangular lattice. For the physi- 
cally more relevant case of quenched disorder, there are 
approximate calcjtiliations using an effective-field theory 
with correlations,!!!! which point to the (expected) weak- 
ening of the tricritical behavior due to the presence of 
disorder. 

In the present work, we first analyze the temperature- 
anisotropy (T x D) phase diagram of the Curie-Weiss 
version (mean-field limit) of the spin Hamiltonian given 
by Eqs. (^) and (||). Depending on the concentration 
p, there appear novel transition lines and multicritical 
points. To include the effects of thermal fluctuations, we 
then resort to a standard self-consistent Bethe-Peierls 
approximation (which is analogous, in the case of the 
corresponding uniform model, to performing an exact 
Bethe-lattice calculation). In section II, we present the 
Curie-Weiss results. The Bethe-Peierls approximation 
is discussed in section III. In the appendices, we report 



an exact solution in one-dimension as well as some low- 
temperature expansions to supplement the numerical re- 
sults from the Bethe-Peierls approximation. 



II. CURIE-WEISS VERSION 

The Curie- Weiss version of the mixed-spin Ising model 
is given by the Hamiltonian 

^ = -^E-^E^. + E^.^^ (4) 

ieA jeB jeB 

where the sums are over all sites belonging to each one 
of the sublattices. 

For a given disorder configuration {Dj}, we calculate 
the partition function by performing a partial trace over 
the set of spin variables {Sj}. In the thermodynamic 
limit, we use the saddle-point method and average over 
disorder to obtain the free energy functional 



1 

2p 



ln2--(l 



a) ln(l 



^ J p{Db) In [1 + 2e-^^« cosh(/3J(T)] dDe. 



(5) 



From the minimization of '^{cr) with respect to cr, we 
have the A-sublatticc magnetization. 



tanh 



(3J / p{Db)- 



2 sinh(/3 Jcr) 



2 cosh(/3 Jo- 



-dDf 



, (6) 



where the random variable Db satisfies the probability 
distribution in Eq. (^. We can now calculate various 
expectation values. For example, we have 



Q = j p{DB){Sl)dDB 

2 cosh(/3 Jcr) 



p{Db)- 



-dDB. 



elSDB + 2 cosh(/3 Jcr, 
The critical line comes from the condition 
92 



= 



cr=0 



1 



(7) 



(8) 



where A = f3D and K = (3J . The thermodynamic sta- 
bility of the critical line depends on the sign of the fourth 
derivative of ^'(cr) at cr = 0. There is then the possibility 
of a tricritical point, given by the additional condition 



= 



3 -I- 9p + ^9 - I86p + I77p2 



cr=0 



(9) 0: 



The stability of the tricritical point is determined by 



> 



■P<Pn 



0.04485 . 



(10) 



which means that the tricritical behavior is suppressed by 
disorder concentrations larger than approximately 4.5%. 

In Fig. 1^, we plot some D x T phase diagrams for a 
set of typical values of the concentration p. In the uni- 
form case {p = 0), there is just a tricritical point, Pj. 
For < p < Prn — 0.04485 . . ., a tricritical point still 
exists (see Fig. |], for p = 0.04). However, at low temper- 
atures and sufficiently large values of D, there appears 
a low-density (Q — > p as T 0) ferromagnetic phase, 
which we call ferro-II phase; at fixed values of D, an 
increase of temperature induces a second-order transi- 
tion from the ferro-II to the paramagnetic phase. This 
transition is represented by a critical line that meets the 
first-order line at a critical end point. Pee- This critical 
end point separates the first-order line into two distinct 
regions: (i) at higher temperatures, there are transitions 
between the usual, high-density (Q — > 1 as T — > 0) fer- 
romagnetic phase (ferro-I) and the paramagnetic phase; 
(ii) at lower temperatures, the transitions are between 
the ferro-I and ferro-II phases, the first-order boundary 
ending at a finite-temperature critical point. Pes- 

For p„, = 0.04485 ...< p < 3/59 = 0.05084 . . ., the 
tricritical point is replaced by a critical end point and a 
simple critical point, separated by a first-order transition 
line between the ferromagnetic phases (see inset in Fig. 
I, for p = 0.05). 



(7=0 



For p > 3/59, the critical line is fully stable (see Fig. 
|l|, for p — 0.08). However, for p < 0.1, there still exists 
a small finite-temperature region where there are (first- 
order) transitions between the ferromagnetic phases. 
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III. BETHE-PEIERLS APPROXIMATION 

To give an estimate of the effects of thermal fluetua- 
tions, which are not accounted for by the Curie- Weiss 
calculations, we now resort to a standard self-consistent 
Bethe-Peierls approximation. As the model is defined on 
a bipartite lattice, we have to consider two distinct clus- 
ters of coordination q (see Fig. ^). In one of them, which 
we call A, the central site is occupied by a spin cr = i, 
connected to q spins of the S = 1 type. In the other clus- 
ter, called B, there is a central S = 1 spin surrounded by 
q spin-i variables. According to the standard prescrip- 
tion of the Bethe-Peierls approximation, we assume that 
the boundary spins of cluster A are under the action of an 
effective magnetic field Hb and an effective crystal field 
D, while the boundary spins of cluster B are in an ef- 
fective magnetic field Ha- The crystal field acting on the 
central site of a cluster _B is a random variable Db- We 



also consider external magnetic fields, Ha and Hb, acting 
on the central sites of clusters A and B, respectively. 

For the sake of simplicity, we assume that the same ef- 
fective crystal field D acts on all boundary spins of clus- 
ter A (and impose self-consistency between both thermal 
and disorder averages associated with the two clusters). 
In a more refined approach, we might introduce different 
effective crystal fields to mimic the extended disorder of 
real materials. However, we expect that, in the absence 
of external magnetic fields, the assumption of a single D 
is reasonable, at least in the paramagnetic phase, where 
there is no long-range order; in particular, we expect that 
the paramagnetic phase boundaries obtained by the two 
approaches are equivalent within the Bethe-Peierls ap- 
proximation. We will see that this is indeed supported 
by the calculations. 

The partition functions associated with the two clus- 
ters are given by 



ZA = e^ 



1 + 2e~^ cosh(7i3 -I- K) 



l-|-2e"'^cosh(7B - K) 



(11) 



and 



Zf 



[2cosh(7A)]'' + e-^^ {e^^ [2cosh(7A + K)f + e^'^^ [2cosh(7A - K)]'^} 



(12) 



where 7 — /3h, A = /3D and K — (3J . The effective 
fields 7^1, 7b and A are determined by the consistency 
equations 



and 



q d^B 

IdlnZA 
q dA 



PyDB)^^^ dDB, 



dAf 



(15) 



where (• • •) and [• • ■]^^ indicate thermal and disorder av- 
erages, respectively. We point out that the introduction 
of the effective crystal field D is essential to achieve con- 
sistency between the equations for the two clusters. 

In order to analyze the critical behavior, it is conve- 
nient to choose the magnetization a, the temperature T, 
and the external fields Hb and Db^ as the independent 
thermodynamic variables. Thus, the external field Ha is 
written as a function of those variables, and the second- 
order transitions in zero external field (Ha — Hb = 0) are 
given by 



djA 
da 



= 



Qo = 



1 



(g - 1)2 tanh^ if 



(16) 



where the derivative is taken for fixed values of the re- 
maining independent variables, and 

2 cosh"? K 



2 cosh^ K 



dDe. (17) 



To calculate the derivative in Eq. (^6|), we take the (im- 
plicit) derivative of the consistency equations with re- 
spect to cr, imposing the condition a = and eliminating 
the derivatives involving S, Q and the effective fields. 
Also, for (T = 0, we have S = ja = Jb — 0, since those 
variables are odd functions of a for Ha = tiB =0. Thus, 
the consistency equation for Q leads to 



= In I 2 

and the final result is 
1 + [2(g 



1 - 



Qo 



■ cosh K 



d^A 
da 



1) - g'l Vo + {q- lyvi 



1 + iq- 2)Vo + {q~ 1)2V? 



(18) 



(19) 



where Vq = Qq tanh^ K. For q = 2, Eq. ( |l9| ) reproduces 
the exact one-dimensional expression for the ^-sublattice 
susceptibility (see Appendix ^). In fact, for g = 2, it is 
not difficult to check that we regain all the exact one- 
dimensional results. 

It is easy to see that, in the uniform case, correspond- 
ing to p{Db) = 5{Db — D), the critical line is given by 



A = In |2 (cosh ^ [q{q - 2) cosh^ K - [q - 1) 



(20) 
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in agreement with the jpsults from the Bethe-latticea and 
the chister-variationalQ calculations. 

Using the binary distribution in Eq. (|^), we obtain 



Qo=P 



2 cosh« K 
1 + 2 cosh« K 



+ il-p)- 



2 cosh« K 



1^ + 2 cosh'' K ' 



Therefore, the critical line is given by 

e- = 2ii^4^cosh^i., 



where 



0W = 



1 cosh-" K 



(21) 



(22) 



2 cosh* K , , 
V . , ......a.. - (23) 



(g - 1)2 cosh^ if - 1 '^l + 2cosh«if 



In the r — !• [K — > oo) limit, we have 



29-1 1 _ 1) 



which has a real solution for A if 
1 - p((? - 1)2 > =^ p < Pc, 



(9-1)^ 



(24) 



(25) 



This last result should be anticipated for a Bethe lat- 
tice, as we can see from the following arguments. Con- 
sider a Cayley tree where the sites belonging to every 
other shell, say those on odd-numbered shells, are oc- 
cupied with probability p, while the remaining sites are 
always occupied. If q is the coordination of the tree, the 
average number of paths from the root (shell 0) to the 
first shell is given by p(g — 1), while we have p{q — 1)^ 
paths from there to the second shell. According to this 
reasoning, the average number of paths from the root to 
the (2n)th shell is given by p^i^q — 1)2". In order to have 
at least one path to the surface of the tree (n — > cx)), it 
is required that -pi^q — 1)^ > 1, which is just the condition 
in Eq. (^5|). This result, together with the reproduc- 
tion of the exact one-dimensional solution, would sug- 
gest that the present treatment also gives exact results on 
the Bethe lattice even in the presence of disorder JJow- 



ever, as remarked in previous similar treatments, EjiIj this 
works for the paramagnetic phase only, because only then 
it is correct to assume that all boundary sites are under 
the action of the same (zero) effective field. The exis- 
tence of a percolating cluster, which we do not take into 
account in this treatment, prevents this approximation 
from still giving correct results for the ordered phases. 

We now consider Eq. (p^, in the infinite coordination 
limit (q ^ oo, if ^ 0, qK = K). We then have 



e^ = 2 



(if - l) - ipif 2 



1 - fpif 2 



(26) 



which agrees with Eq. 
the model. 



for the Curie- Weiss version of 



The tricritical points are determined by Eq. (^^ sup- 
plemented by the condition 



= 0, 



a=0 



which is equivalent to 

2g2 -lOq + 6 
{q- l)5tanh2is: ^ 

where Wq is given by 



3qWo tanh^ K 



(g-2)(g-3) 

{q-ir ■ 



Wo = / piDe) 



2 cosh'' K 



+ 2 cosh'' K 
The tricritical points are stable if 



dDf 



(27) 



(28) 



> 0. 



To calculate this derivative, we again take implicit deriva- 
tives of the consistency equations (up to fifth order) with 
respect to cr, at cr = 0, and eliminate all derivatives in- 
volving S, Q and the effective fields. Unlike the previous 
analysis, we have not been able to obtain closed- form 
expressions for the stability condition of the tricritical 
point, but it is not difficult to perform a number of nu- 
merical calculations. 

For the uniform model, we have Wq — Qq. Therefore, 
Eq. (p^) takes the form 



tanh K = 



1 



5q~3 
1^' 



(29) 



which is again identical to the result obtaijaed from 
the Bethe-latticeH and the cluster-variationalB calcula- 
tions. Notice that this equation has real solutions only if 
q > 4.561553 . . .. Thus, the Bethe-Peierls approximation 
does not predict a tricritical point for the square lattice 

(q = 4). 

For the binary distribution in Eq. (|^), we have 



1-p 



1 2 cosh'' K 
1 + 2 cosh'' K 



(30) 



In the infinite-coordination limit we can write 

2' 



if4 



1 



P 



1-p 



which leads to the equation 
p 1 . 4 ~, 



1 



1-p 



1 



2 ~, 
1 - -if 2 
3 



lif^ 
9 



(31) 



0, (32) 
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at the tricritical point. Indeed, one of the solutions of this 
equation corresponds to Eq. (g), for the Curie- Weiss ver- 
sion of the model, while the other solution represents a 
thermodynamically unstable situation. 

In Table ^, for various values of the coordination num- 
ber q, and using the binary distribution given by Eq. 

, we give the corresponding values of the concentra- 
tion pm, at which the tricritical point becomes unstable, 
and the critical percolation concentration p^r- We see 
that, for q < 10, the tricritical behavior is suppressed 
for pm < Per, while, for q > 11, that suppression occurs 
for pm > Per- As shown in Table |[ Pm increases with q, 
which indicates that disorder is more effective for small 
coordination numbers. 

As the effects of binary disorder strongly depend on 
the coordination, we now discuss the phase diagrams for 
the typical cases. 

For g = 3 and 4, there are no tricritical points. The 
D X T phase diagram displays just a fully stable critical 
line. The main effect of disorder is to make the param- 
agnetic phase unstable at T = 0, regardless of the value 
of D, for p larger than the critical percolation concen- 
tration Per- The phase diagrams in Fig. |[ for g = 3, 
are in qualitative agreement with the exact results for 
the honeycomb lattice (which is also three-coordinated) 
under annealed disorder .IMj At T = 0, there is even quan- 
titative agreement with the value of the critical crystal 
field at pcr, given by Dcr = 5J/3, although of course 
this agreement does not extend to the value of pcr it- 
self. Our results for q = 3 and q = 4 are also in qual- 
itative agreement with those obtained by a real-space 
renormalization-group approach for the two-dimensional 
Blumc-Emery-Grifhths model in a random crystal field.Eil 
For 5 < q < 10, the concentration pm above which 
the tricritical point becomes unstable is lower than pcr- 
For p < p„i, disorder dcpressses the tricritical temper- 
ature, and shortens the first-order transition line. For 
Pm < p < Per, the tricritical point is replaced by a crit- 
ical end point. Pee, and a simple critical point. Pes, as 
in the Curie-Weiss version of the model. However, the 
paramagnetic phase is stable at T = if £> > qJ, and the 
first-order line reaches D — qj a.t T — 0. As p increases, 
first the critical end point Pee and then the simple critical 
point Pes reach the T = axis, at values of p which can 
be determined by a low-temperature expansion of the 
consistency equations (see Appendix P). In Fig. |[ we 
plot the D X T phase diagram for q = 6 and p — 0.011. 
To determine the first-order lines shown in that figure, 
we numerically solve the consistency equations to obtain 
the conditions Ha (ci) = Ha {<J2) = and 

J hA (cr) da = 0, (33) 

which correspond to a Maxwell construction. 

For q > 11, we have pm > Per, so the behavior of the 
system is quite similar to the predictions of the Curie- 



Weiss version of the model. 



IV. CONCLUSIONS 

We performed detailed calculations for the phase di- 
agram of a random-anisotropy mixed-spin Ising model 
both in the mean-field limit (Curie-Weiss version of the 
model), in which thermal fluctuations are neglected, and 
according to a standard self-consistent Bethe-Peierls ap- 
proximation (which turns out to be exact in one dimen- 
sion). For a binary distribution of crystal fields, we ob- 
tained closed-form expressions for the critical lines and 
the location of the tricritical points. Depending on the 
concentration p, the mean-field results for the D x T 
phase diagrams predict novel first-order lines and mul- 
ticritical points (besides a ferromagnetic region, for all 
values of the crystal field, extending down to the lowest 
temperatures). The Bethe-Peierls approximation shows 
that this additional ferromagnetic region is suppressed 
for concentrations below a certain percolation threshold. 
Also, the Bethe-Peierls results point out to the absence of 
a tricritical behavior for lattices with coordination g < 4. 
All results reported in this paper are in agreement with 
general predictions for the effects of disorder on first- 
order transitions and multicpitical points (for a recent 
review, see a paper by CardyllJ). 
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APPENDIX A: EXACT SOLUTION IN 
ONE-DIMENSION 

For an open chain with N + I sites {N even) , and in 
zero external field, the Hamiltonian of the mixed-spin 
Ising model can be written as 

N/2 N/2 

i=i i=i 

Given a disorder configuration {D} = {-Di, . . . ,-DAr/2}, 
we perform a partial trace over the spin variables {Sj} 
to write 

{•^y {s} 

N/2 

= E n {1 + cosh[K{a, + a,+i)]} , (A2) 

where K — (3, J and Aj = PDj. Introducing a prcfactor 



(D 



5 



^2 = (1 + 2e-^^) [1 + 2e-^^ cosh(2if)] (A3) After some algebraic manipulations, for j < k, we have 



and an effective interaction X,-, such that 



1 + 2e~ 



cosh(2iv:) 



1 + 2e-^j 



(A4) 



the partition function can be written as the factorized 
form 



N/2 



{a} J = l 
N/2 

= J]^ 2 [l + 2e-~^^ cosh^ K] . 



(A5) 



From Eq. (A5), we obtain the thermal average 
dlnZ 2e~^^ cosh^ K 



{D} 



9A, 



1 + 2e-'^^ cosh^ K ' 



(A6) 



which depends on the value of the crystal field on the 
jth site only. Since we are considering a nearest-neighbor 
one-dimensional model in zero field, the thermal averages 
{Sj) and (aj) are zero. Performing the disorder average, 
we obtain the expectation value 

r ""^ r 

i=l 

(A7) 

For a given disorder configuration, the magnetic sus- 
ceptibilities of the a and S sublattices are given by 

XAD} = I jim J2 E (A8) 



T N^oo N + 2 



j=l k=l 



and 



N/2 N/2 



j = l k=l 

The two-spin correlation functions, 
^ ^ W) {S} 

and 

-!3H 



(A9) 



(AlO) 



{S. 



^^\o} = ^^llT.S^S,e-P". (All) 
^ ^ {S} 

can be calculated if we introduce the transformation 

Tj = UjUj+i with To = (Ti. (A-12) 



^ ^ 2smh^K 



(-.^^■){i.} = ll,A.+2cosh^if 



(A13) 



and 



n 

\i=hk 
fe-1 

n 



sinh 2if 



s^' + 2 cosh^ if 



2 sinh^ K 
3^- 2 cosh^ K ' 



(A14) 



from which we obtain the expectation values 

„ N/2 

9a{\k~j\)^ I {ajak)^jj^'Wp{Di)dD^ 



and 



(gtanh^ Kf~'^ 



N/2 



(A15) 



= Q{Q tanh^ is:^ 



(A16) 



which depend on the distance between sites j and k. The 
expectation values of the susceptibilities are given by 



1 

T 



1 + 2^3, (r) 



and 



1 

T 



Q + 2j29s {r) 



1 1 + Q tanh^ K 
Tl-Q tanh^ X 

(A17) 

Ql + Q tanh^ 
Tl-Q tanh^ X ' 

(A18) 



where Q is determined by Eq. (A7) 



APPENDIX B: LOW-TEMPERATURE 
EXPANSION 

For the binary distribution, in the low-temperature 
limit {K = (3J ^ 1), if we neglect terms of order 
exp {—2K) and higher, the consistency equations (p^- 
( p^ ) for cluster A lead to the expressions 



7/1 = — hi — hi -— , 

' 2 1 - (T 2 1 - C_ ' 



5 = 



(Bl) 



(B2) 
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and 



where 



c±^ 



1 -cr 



= ±7B 



(B3) 
(B4) 



For cluster B, we have 

a = ptanh {qja) 

t{-ja) tanh(7A) + S tanh(g7^) 



+ {l-p)- 



S ^ p tanh (qja) + (1 - p) 
Q^P+il-p) 



5 



(B5) 



t(7a) + S 
S 



where 



and 



t{x) = 



TilA) - 

5 = exp{qK - A), 
2" 



(1 + tanh a;)' + (1 - tanhx)'' 



tanh (7^) , (B6) 
(B7) 

(B8) 
(B9) 



Solving Eqs. ( p32D and (|B3|) for C'± in terms of cr, S 
and Q, and using Eqs. (B5)-(|B7|), we can write Eq. (Bl) 
in the form 



i^lni±^ 

2 1 - (T 



(BIO) 



where 7a (f) is determined from t he so lution of E q. (pa). 
Notice that, according to Eqs. ( B10| ) and (B5), 7a(ct) 
and 7^(0") depend on the temperature through the pa- 
rameter 6 only. As r — !■ 0, this parameter goes to zero (if 
D > qj) , or infinity (if D < qj) , except in the vicinity of 
the point Po with coordinates D = qJ, and T = 0, where 
S can assume any value. 

Since the equation of state ( B10| ) becomes asymptot- 
ically exact as r — > 0, it can be used to determine the 
values of p at which the critical end point and the simple 
critical point reach Pq, and thus disappear. To do that 
calculation, we impose the conditions 



7A(cre) 



djA 
da 



d^lA 



0, 



(Bll) 



from which we obtain the values of a,, , 5^ and Pe at which 
the critical end point reaches Pqj and the conditions 



djA 
da 



-1A{a)da = {), (B12) 



which give the corresponding values as , Ss and Ps for the 
simple critical point. 
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FIG. 1. Phase diagrams of the Curie-Weiss version for typ- 
ical values of the disorder concentration p. 
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FIG. 2. Clusters used in the Bethe-Peierls approximation. 
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FIG. 3. Phase diagrams for coordination q = 3 according 
to the Bethe-Peierls approximation. 
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FIG. 4. Phase diagram for coordination q = Q and disor- 
der concentration p = 0.011 according to the Bethe-Peierls 
approximation. 

TABLE I. Values of the critical percolation concentration 
Per and the concentration pm at which the tricritical point 
becomes unstable, as functions of the coordination q in the 
Bethe-Peierls approximation. 



q Per Pm 

~5 6.25 X 10"'' 7.4161 x 10"* 

6 4 X 10"^ 2.0454 x 10"^ 

10 1.2346 X 10"^ 9.8265 x 10"^ 

11 1 X 10"^ 1.1665 x IQ-^ 
20 2.7701 X 10"^ 2.3001 x 10"^ 
100 1.0203 X 10"* 3.9707 x 10"^ 
oo 4.4850 X 10"^ 
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